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ABSTRACT 


A  microscopic  model  for  the  metal-solution  interface  is  presented t 
which,  for  the  first  time,  includes  both  a  non-ideal  treatment  of  the 
metal  and  a  molecular  model  for  the  solution  side.  The  metal  is  described 
by  a  jellium  including  electron-ion  pseudopotential.  The  solution  is  treated 
as  a  mixture  of  hard  spheres  with  point  dipoles  (solvent  molecules)  or  point 
charges  (ions).  The  statistical  mechanics  is  solved  using  the  mean  sphe¬ 
rical  approximation.  No  a  priori  separation  is  introduced  between  compact 
and  diffuse  layers.  A  simplified  version  of  the  model  is  applied  to  the 
case  of  Hg-DMSO  systems. 


INTRODUCTION 


The  classical  theories  of  the  ideal  polarized  electrode  are  based  on 
two  premises  : 

a)  The  metal  has  no  direct  effect  on  the  differential  capacitance  C» 
although  it.  does  determine  the  potential  of  zero  charge  and  affects  the  orien¬ 
tation  of  solvent  molecules.  From  an  electrostatic  point  of  view,  every  metal 
acts  just  the  same,  the  only  interaction  being  that  of  ideal  images. 

b)  The  solution  side  has  two  layers  :  the  external,  or  diffuse  layer, 
which  can  be  described  by  such  simple  models  as  the  Gouy— Chapman  (G.  C.) 
theory,  and  the  internal;  or  compact  layer,  which  in  the  absence  of  specific 
adsorption  is  formed  only  by  solvent  molecules. 

These  two  premises  are  now  being  challenged.  The  important  progress 
in  the  theory  of  metals  and  metal  surfaces  in  recent  times  makes  it  hard  to 
admit  that  the  metal  is  totally  unresponsive  to  changes  in  the  electrolyte. 

The  surface  charge  distribution  in  the  metal  must  respond  to  local  changes  in 
the  electrostatic  fie^d  produced  by  modifications  in  the  ionic  and  solvent 
distributions.  Surely  this  will  produce  contributions  to  the  differential 
capacitance,  that  will  be  different  for  each  metal,  since  the  response  to  an 
external  field  is  different  for  each  metal,  and  depends  on  details  of  the 
electronic  structure  near  the  surface.  In  previous  work  [1]-  [3]  an  estimate 
of  these  effects  was  given.  The  solution  side  was  represented  by  a  laye£  of 
dielectric  material.  In  a  recent  publication,  W.  SCHMICKLER  [4]  has  discussed 
a  similar  model.  His  conclusions  about  the  role  of  the  metal  are  similar  to 
ours. 

On  the  solution  side,  recent  advances  in  the  treatment  of  molecular 
solutions  [5] [6],  that  is  solutions  in  which  the  solvent  is  not  treated  as  a 
continuum,  have  made  it  possible  to  treat  the  statistical  mechanics  of  the 
solution  side  of  the  electrode  interface.  Although  the  approach,  in  principle, 
can  handle  realistic  models  of  solvents  and  ions,  the  initial  work  has  dealt 
with  a  model  in  which  the  solvent  is  represented  by  hard  soheres  with  point 


dipoles  [7].  This  work  was  revealed  two  important  deficiencies  of  the  classic* 
treatment  : 

a)  The  separation  of  the  inner  and  outer  layers  for  the  solvent 
orientation  process  (or  better  polarization)  implies  a  violation  of  Maxwell's 
equations.  The  dipoles  of  the  solvent  are  oriented  wherever  there  is  an 
electric  field,  that  is,  throughout  the  entire  electric  double  layer. 

local 

b)  The  use  ofVdielectnc  constants  is  also  unwarranted.  Although 
this  is  a  much  more  subtle  problem,  recent  theoretical  work  [8]  has  shown 
that  dielectric  constants  have  a  clear  meaning  only  for  large  systems  (ther¬ 
modynamic  limit)  and  not  at  the  molecular  level.  The  polarization  is  the  local 
variable  that  is  required. 

In  this  work  we  present  a  comprehensive  theory  in  which  all  of  these 
concepts  are  included.  In  section  2  we  review  the  models  of  the  metal  and 
solution  sides  of  the  interface.  The- potential  drop  and  capacitance  are  then 
computed  using  a  simple  GMSA  based  model  which  has  the  obvious  advantage  of 
analiticity  and  simplicity. 

In  section  3  we  discuss  the  electronic  density  profile  of  this  model. 
We  consider  then  a  model  in  which  the  physical  parameters  are  those  of  an 
alkali-ion  in  DMSO  in  the  neighborhood  of  a  mercury  electrode. 

In  our  model  the  point  of  zero  charge  is  not  the  point  of  zero 
potential  (which  is  the  case  in  the  restricted  primitive  model)  :  we  discuss 
also  the  change  in  the  surface  potential  of  the  metal  induced  by  the  presence 
of  the  solution.  The  influence  of  some  of  the  parameters,  and  detailed  compa¬ 
rison  between  Hg  and  Ga  electrodes  will  be  discussed  in  future  work. 


In  the  last  section  we  examine  possible  improvements  in  our  work- 


2.  MODEL  FOR  THE  SURFACE  INTERACTIONS 


Consider  a  plane,  ideally  smooth  electrode,  z  is  the  distance 'from 
the  metal  surface.  The  bulk  metal  is  then  given  when  z  -*■  -  «  .  For  the 

metal  we  use  a  jellium  model  (with  inclusion  of  the  electron-ion  pseudopoten¬ 
tial).  The  details  are  given  in  [9],  but  for  the  sake  of  completeness  we  will 
give  an  outline  of  this  work.  The  metal  ions  are  taken  as  a  continuum  of 
density  p  for  z  <  0  .  In  other  words 

p(z)  -  p  0(-z) 

where  6(x)  is  the  Heaviside  function 

0(x)  *1  x  >  0 

0(x)  -  0  x  <  0 

The  electronic  gas  is  free  to  move,  but  its  equilibrium  distribu¬ 
tion  n(z)  will  be  the  one  that  minimizes  the  total  surface  energy  of  the 
system.  The  electron-ionic  background  interaction  is  given  by  a  pseudopotential 
The  properties  of  the  electron  gas,  such  as  the  kinetic,  exchange  and  correla¬ 
tion  energies  are  computed  using  the  local  approximation. 

On  the  solution  side  the  ions  are  represented  by  hard  spheres  of 
diameter  a.  and  charge  ±  e::,  and  bulk  density  p..  The  solvent  is  also  repre¬ 
sented  by  hard  spheres  of  diameter  cr^,  point  dipole  u  ,  and  density  p^. 

We  assume  that  no  chemical  reaction  can  take  place,  no  miscibility 
between  the  phases  [10].  However,  this  model  does  not  exclude  the  existence 
of  a  surface  layer  in  which  all  the  particles  are  present.  In  fact,  the 

existence  of  this  layer  is  cons  is  cant  with  the  contemporary  theories 
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of  the  metallic  surfaces.  Any  model  which  does  not  include  this  layer  leads 
to  predictions  which  are  incompatible  with  experiments  [1][11]. 

In  our  model,  the  electrons  can  diffuse  into  the  ionic  solution. 
In  actual  practice,  calculations  show  that  both  in  the  interface  with  the 
vacuuL12,1  or  with  a  dielectric  film^,  the  electron  density  spill  over  is 
very  small.  However  the  change  in  the  spill  over  with  the  electrode  charge 
will  give  a  direct  contribution  of  the  metal  to  the  differential  capacitance 

[1][4]. 

The  metal  electrons  in  the  solution  are  subject  to  electrostatic 
interactions  (electron-ion  and  electron-dipole)  but  also  a  repulsion  due  to 
the  core  electrons  of  ions  and  solvent  molecules.  These  are  represented  by 
the  Harrison  repulsive  potential  [9]. 

There  are  several  parameters  in  our  model  that  are  fixed  a  priori, 
but  which  should  come  out  of  a  first  principles  calculation.  One  of  these 
quantities  is  the  distance  of  closest  approach  of  the  ions  and  solvent  to 
the  plane  z  ■  0  (the  metal  surface).  We  will  assume  that  it  is  just  the  sum 
of  the  molecular  radius  of  the  solvent  or  ions  (which  for  this  simple  model 
are  considered  to  be  equal)  and  the  radius  of  the  metallic  ion.  In  reality, 
this  parameter  should  arise  from  a  self-consistent  density  functional  calcu¬ 
lation  [13  ] . 

In  short,  we  will  assume  that  there  is  an  ideally  smooth  wall, 
located  at  z  ■ 6  ,  which  cannot  be  penetrated  by  either  ions  or  solvent  mole¬ 
cules.  The  distances  of  closest  approach  are 
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solvent  molecules 

There  is  also  the  problem  of  surface  smoothness  and  structure. 

From  recent  experimental  evidence,  it  is  known  that  certain  metal  surfaces 
are  subject  to  structural  changes  during  chemisorption  processes  -  In 
these  cases  the  surface  overlayer  changes  its  structure.  It  is  quite  concei¬ 
vable  that  in  an  intrinsically  more  labile  environment  such  as  that  of  liquid 
metals,  the  metal  overlayer  would  also  respond  structurally  to  changes  in 
potential . 


Furthermore,  the  question  of  surface  ripples  (capillary  waves)  at 
the  interface,  and  its  effects  on  the  structure  and  thermodynamics  of  the 
electrode  is  not  discussed.  We  hope  to  address  these  and  other  questions  in 
future  work. 

Summarizing,  the  ions  and  dipoles  in  the  solution  do  not  see  simply 
a  hard,  charged  surface,  but  rather  the  sum  of  a  charged  surface  (located  at 
6  +  a/2  )  and  a  distribution  of  electron  density  beyond  the  plane  of  closest 
approach. 


Even  this  simple  model  is  not  easily  amenable  to  a  complete  nume¬ 
rical  calculation.  We  will,  therefore,  in  this  first  communication  ,  use  a 
simplified  theory  that  has  the  advantage  of  giving  explicit  analytic 
expressions  for  the  physically  relevant  quantities. 


3.  A  SIMPLIFIED  MODEL 


For  simplicity,  we  shall  discuss  only  the  restricted  model  in  which 
a^  ■  a+  ■  a_  *  a  (Although  not  in  published  form,  the  results  for 

a+  ■  a_  /  (Jj  are  also  available).  We  neglect  also  the  effects  of  the  electron 
density  on  molecules  and  ions  in  the  solution.  We  shall  leave  the  discussion 
of  the  validity  of  this  assumption  for  the  future1:  The  solution  is  represented 
by  the  MSA  [5] [6]  which  is  only  qualitatively  correct.  Near  a  charged 
interphase  a  simplified  theory,  valid  presumably  for  dilute  solutions,  has 

been  developed  by  Blum  and  Henderson  [7],  In  this  theory,  the  dipolar  interac- 

% 

tions  are  decoupled  from  the  ionic  distributions,  which  are  given  by  the 
modified  Gouy-Chapman  (MGC)  theory.  We  are  thus  satisfying  the  electrostatic 
part  of  the  contact  density  theorem.  The  simplified  model  of  the  interface 
is  given  on  Fig.  !. 


The  charge  density  profile  is  given  by 
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where  <  is  Debye's  inverse  length  and 


dj  »  6  +  a/2 


The  total  excess  charge  of  the  solution  side  is 

ao 

%  -  /  dz  -  -  qM  (?) 

di 

The  electronic  density  profile  is  postulated  to  be  of 
form  [1][9]  : 


n(z)  -  n  [I  -  -i  ea^z”*o^] 
-  (n/2) 


z  <  z 


where  n  is  the  bulk  metallic  electron  density.  .Clearly,  the  total  excess 
charge  on  the  metal  side  is 


qM  *  /  dz  Cp© (— z>  -  n(z)]  -  - 


This  last  relation  also  defines  z  . 

o 


For  a  given  value  of  a  ,  and  q.,,  the  charge  and  polarization  (P(z)) 
profiles  are  computed  from  the  theory  [7].  These  in  turn  yield  the  potential 
V(z)  and  electric  field  E(z),  which  are  obtained  by  solving  Poisson's 


equation 
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The  potential  V(z)  is  continuous  in  all  space.  Because  of  the 
discontinuities  in  the  density  profiles,  V(z)  has  a  different  form  in  the 


regions 


I)  -  08  <  z  <  0 

II)  0  <  z  <  dj 

III)  dj  <  z  <  oo 
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The  soluCion  of  (5)  is  straightforward,  and  is  explained  in 
Appendix  I.  The  total  potential  drop  across  the  interface  is 
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The  integral  in  this  expression  is  explicitly  known.  The  reader 
is  referred  to  Ref.  7  (eq.15)  for  the  derivation.  The  result  is 
» 
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The  Laplace  transfers  P(s)  is 


'U 

P(s) 


4ir  / 
d 


dz  e  SZ  P(z) 


1 


47ryp ,  By4frq 


( 


M 


I 


rr  rr 


(<  +  s)  A(s) 


(7) 


b: 

{. 


b„ 

<o  3 


-t  rr  *  2 


12  B 


Here  is  the  bulk  number  density  of  solvent  molecules.  The 
coefficients  83,  8g  and  B12  are  functions  of  the  bulk  dielectric  constant  e 
of  the  solvent,  and  are  given  through  an  auxiliary  parameter  A 
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The  auxiliary  parameter  X  is  computed  from  the  equation 

X2(!4-X)4 
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The  functions  $,(s)  and  A(s)  are  defined  by 


chj  Cs)  =•  l/(Sa)2  [l  -  SC  -  e"sa] 
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Substituting  into  equation  (6)  yields 
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at  the  point  of  zero  charge  (zq  ■  0) ,  this  expression  simplifies  to 


The  potential  drop  across  the  interface  is  a  function  of  only  n 

and  a  ,  which  are  metal  parameters.  This  however  does  not  mean  that  the 

potential  AV  is  equal  to  the  surface  potential  of  the  metal-vacuum  interface 

X^t  since  a  is  also  a  function  of  the  solution  parameters.  In  the  limit  a-** 

the  electronic  profile  tends  to  n(z)  ■  n  9(z-zq),  and  then  Eq.(Il)  is  not 

the  same  as  the  expression  given  in[14].  There  are  two  reasons  for  this 

2 

discrepancy.  First,  the  term  2imzo  does  not  exist  in[14],  because  now  the 
charge  is  not  localized  at  the  plane  z  -  0.  Secondly,  the  third  term  in  (11) 
is  not  that  of[14].  However,  if  dj  ■  o/2  then  they  are  identical.  This 

means  that  an  additional  capacitance  term  appears  because  of  the  shift  5 
in  the  metal  background  profile. 

% 

The  differential  capacitance  of  the  system  is 
1  3AV  3AV  3AV 


using  (11),  we  obtain 
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In  Chis  expression  Che  first  term  of  the  right  hand  side  is  Che 
result  of  the  MSA  [7] [8].  If  the  solvent  molecules  are  shrunk  to  point 
dipoles  <(7^  -*■  0),  one  should  recover,  from  (14)  the  differential  capacitance 
of  the  modified  Gouy  Chapman  theory.  Otherwise  £14]  one  can  define  the  effec¬ 
tive  dielectric  constant. 
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which  then  yields  the  MGC  result  with  the  effective  dielectric  constant  £•  We 
must  remark  that  in  spite  of  the  formal  res semblance,  the  physics  is  diffe¬ 
rent,  since  the  solvent  is  polarized  throughout  the  entire  electric  double 
layer. 


The  second  term  of  (14)  arises  from  the  fact  that  the  electrode 
charge  is  not  localized  at  z  ■  0.  The  next  term  is  due  to  changes  in  the 
dipolar  surface  layer  of  the  metal  as  a  function  of  the  excess  charge.  And 
the  last  term  reflects  the  fact  that  the  distance  dj  is  also  a  function  of 
the  excess  charge.  A  detailed  discussion  of  this  fact  can.be  found  in  [3]. 
At  the  point  of  zero  charge  this  term  vanishes. 


Because  of  the  intrinsic  limitations  of  the  MSA,  we  will  limit 
ourselves  to  small  surface  charge  densities  around  the  point  of  zero  charge. 

Thus,  for  a  given  value  of  dj,  we  only  need  to  compute  a  for  the 
differential  capacitance  C. 
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4.  ELECTRONIC  DENSITY  PROFILE 


The  nns.itz  (3)  or.  the  electronic  density  profile  is  particularly 
convenient  since  n(z)  only  depends  on  the  parameter  a  .  We  will  compute 
this  parameter  by  minimizing  the  surface  excess  energy  Ug .  The  kinetic  energy 
and  the  non-coulombic  part  of  the  electron  -  neutralizing  background  are  not 
affected  by  the  presence  of  the  solution,  and  therefore  we  use  the  expres¬ 
sions  given  in  Ref. [9]  (Eq.  14  with  a  ■  S) . 


tt_ .  +  U  ♦  U.  .  -  -  1.6423  n2/3  +  0.25037  n1/3  +  .006584] 

Km  exc  inn  a 


(16) 


on(log2)/72 


The  pseudopotential  contribution  U  to  the  energy  is  given  by 


(Eq.  10  of  [1]  with  a  •  $)  : 
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where  Aq  and  are  the  pseudo  potential  parameters  of  Heine  and  Animalu 
(see  also  [1]  ). 

The  electrostatic  energy  is  computed  from  the  charge  and  polariza¬ 
tion  profiles,  and  the  potential  V(z)  which  has  been  calculated  in  Appendix  I. 
We  find 
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es 
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The  last  term  consists  of  the  remaining  contributions  to  the 
energy  which  are  not  functions  of  a,  and  thus  irrelevant  to  our  calculation. 

The  term  before  the  last  in  (18)  represents  the  interaction  of  the 
electrons  with  the  local  polarization  of  the  solvent  molecules  in  the  solution. 
This  contribution  is  easily  calculated  from  the  Laplace  transform  P(s)  (Eq.7) 
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There  is  also  another  contribution  due  to  the  repulsive  interactions 
of  the  electrons  and  the  electronic  clouds  of  the  ions  and  solvent  molecules 
[9].  These  interactions  are  represented  by  Harrison's  pseudopotentials*  with 
parameters  A.  (for  ions),  A^  (for  solvent  ).  This  interaction  is  independent 
from  the  ionic  charge  or  solvent  dipolar  orientation,  and  because  of  the  us* 
of  the  MSA,  it  will  depend  only  on  the  total  density  profile  pQ(z).  Since 
all  the  particles  are  hard  spheres  of.  equal  diameter,  a  single  parameter  X 
(a  weighted  average  of  A^  and  A^)  will  be  enough  for  our  purposes.  Ve  get 

UQ  ■  X  /  dr  /  dz'  n(z)  6(r  -  r')  pQ(z')  (20 

-  Vd  +  2Xip* 

pd+2p+ 

where  is  the  dipole  concentration  and  p+  is  the  (+)  ion  concentration. 

The  integration  is  easily  performed 

UB  **  X  /  dz  P0(z)  '  (2D 

dl 

-  ^  eaz°  pQ(a)  (21*) 

*\j 

where  pQ(a)  is  the  Laplace  transform  of  the  total  density  distribution 
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The  functions  pj(a)  and  are  defined  by  Eqs  (9). 
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The  total  surface  energy  is  then 


U  -  IL  .  +  U  ♦  U.  .♦  U  +  U  +  U„ 
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Using  Eqs  (16-22)  we  get.  finally 
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The  coefficients  C^,  C^,  E.  D,  are  given  in  the  Appendix  2. 

In  Equation  (23),  all  the  electrostatic  interaction  terms  with 
the  solution  cancel  at  zq  ■  0,  the  point  of  zero  charge ✓  The  shift  in  the 
surface  potential  of  the  metal,  when  it  is  taken  from  the  vacuum  to  the 
solution  is  due  to  the  electron— solution  repulsive  interactions. 


The  minimization  of  U  yields  then  a  .  This  parameter  is  them 

s 

used  to  compute  the  properties  of  the  electrode  with  the  aid  of  (11) 
and  (19). 


-  i  J  - 


5.  RESULTS 

Clearly  the  MSA  as  a  weak  coupling  parameter  theory  is  not  suited 

to  describe  the  behaviour  of  hydrogen  bonded  solvents,  such  as  water.  However, 

an  aprotic  solvent  like  DMSO  can  be  reasonably  described  by  the  MSA.  Indeed, 

3 

with  the  experimental  values  for  the  density  p.  (14.05  M/dcm  at  room  tempe- 

a 

rature)  and  the  permanent  dipole  moment  u  (  *\*  4  Debye) ,  the  MSA  yields 
e  ■  58  instead  of  the  experimental  value  47.15  (see  for  instance  [8]  for 
the  MSA  expression  of  e  ).  From  (8')  this  gives  for  X  respectively  2.49  and 
2.39. 

From  (14)  we  see  that  the  molecular  contribution  to  the  different 

tial  capacitance  C  depends  only  on  e  ,  while  from  (23)  a  depends  on  X  (or  c), 

Pd  and  u  •  We  have  checked  that  changing  X  does  not  affect  significatively 

the  value  of  a  which  minimized  U  ,  and  hence  the  metal  contribution  to  the 

s 

capacitance.  Therefore  we  have  used  X  -  2.39  (i.e.  £  -  47.15)  in  the  nume¬ 
rical  computations. 

We  now  choose  the  effective  diameter  of  the  DMSO  molecule  so  that 
the  parameter 


n  “  6  pd 


0.45 


which  is  the  value  of  the  volume  fraction  characteristic  of  many  liquids. 

0 

This  yields  a  «  8.81  a.u.  (4.66  A),  which  is  a  little  small,  but  not 
unreasonable  for  this  molecule.  The  ionic  diameters  are  taken  also  to  be 
8.31  a.u.  Fig.  2  shows  the  corresponding  MSA  density  profile  for  thejnole— 
cules  which  have  their  dipole  moment  normal  to  the  interface. 


The  electrode  is  liquid  mercury,  for  which  the  parameters  n,  Rjj 
and  A  are  taken  from  [9].  For  the  mercury-vacuum  interface  we  found 

=  0.901  u.a.  and  a  surface  potential  =  5.34  eV.  Since  we  will  also 
include  some  results  for  gallium,  we  refer  the  reader  to  our  previous  work 
[9]. 

Two  of  the  parameters  of  the  problem  are  free,  and  have  to  be 
selected  using  physical  intuition.  They  are  5  and  X  .  It  seems  realistic 


eo  choose  Che  value  o£  6  between  0  and  R  Che  crystallographic  radius  of 
the  ions  in  the  metal  [9]. 

For  water  the  parameter  X  has  been  estimated  to  be  'V*  15.  For  SMSO 
we  might  expect  an  even  higher  value  since  there  are  54  electrons  in  DMSO 
against  10  in  water.  In  table  1  we  show  the  influence  of  X  on  a  assuming 
that  dj  ■  a/2  or  d^  -  Rc  +  a/2  .  As  expected  the  repulsion  term  with  the 
solution  causes  an  increase  in  the  values  of  o  .  This  is  similar  to  what  ve 
found  previously  [9],  In  fig.  (3)  we  have  plotted  Che  shift  in  the  surface 
potential  of  the  metal  5^  ■  4to  (—j  -  ~)  calculated  at  the  point  of  zero 
charge  as  a  function  °  ao  of  the  distance  dj.  The  corres¬ 

ponding  case  of  gallium  for  X  ■  30  is  also  shown. 

The  repulsion  from  the  solution  is  not  enough  to  keep  the  electrons 
in  the  region  z  <  dj.  There  is  always  a  small  charge  spillover,  -  6q,  in 
the  region  z  >  dj. 


n  -a(z-z) 

IE  e  ° 


In  fig.  (4)  we  show  5q<dj)  for  both  Hg  and  Ga.  The  spillover  fiq 
is  larger  for  Ga  chan  for  Hg. 

In  the  simple  model  presented  here,  this  effect  has  been  neglected 
in  calculating  the  charge  distribution  and  the  polarization  in  the  solution. 
It  is  clear  that  this  approximation  is  less  justified  for  Ga  than  for’Hg. 
However,  the  mathematical  tools  to  include  it  are  available  and  we  will 
indeed  discuss  these  effects  in  a  near  future. 

The  fact  that  there  are  electrons  on  the  solution  side  even  at  the 
point  of  zero  charge,  indicates  that  the  solvent  molecules  must  be  polarized 
at  the  interface  even  at  the  point  of  zero  charge.  Note  that  the  ratio 
5q(Ga)  /  6q(Hg)  is  not  simply  the  ratio  of  the  electronic  densities. 

The  differential  capacitances  for  the  mercury  electrode  are  presen¬ 
ted  in  table  II.  They  were  computed  from  the  curve  AV  -  f(q„)  bv  numerical 
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differentiation.  With  the  parameters  of  this  article  the  zero  charge  capaci¬ 
tance  corresponding  to  an  ideal  metal  would  be 


ci  "  ^  — to - £zr~  '  7-94  vF,ca2 

**  1  ♦  ~  a  + 


distance  d 


We  see  that  the  total  capacitance  crucialy  depends  on  the 
1* 


When  dj  is  large  the  value  of  C  is  greater  than  but  the  effect 
of  the  non-ideality  of  the  metal  remains  small.  On  the  other  hand,  when  the 
distance  of  closest  approach  decreases  (for  instance  if  there  is  interpene¬ 
tration  between  the  metallic  ions  and  the  particles  of  the  solution)  this 

effect  becomes  quite  important  and  the  influence  of  X  is  not  negligible  (for 

2 

comparison  the  experimental  value  if  18  yF/cm  for  LiClO^  ions  [15]  ). 

First  calculations  show  chat  the  differential  capacitance  for  Ga 
would  be  greater  than  for  Hg,  as  observed  experimentaly. 


6.  CONCLUSIONS 


We  have  presented  a  realistic  model  o£  the  metal  solution  interface, 
in  which  the  metal  is  not  an  ideal  smooth  metal,  but  a  jellium  model.  A 
simple  discussion  of  this  model  is  given.  The  results  confirm  earlier  calcu¬ 
lations  [l]-[4]  that  indicate  that  the  metal  makes  an  important  contribution 
to  the  differential  capacitance  C. 

The  solution  side  of  our  model  is  a  mixture  of  hard  spheres  with 
charges  and  point  dipoles.  No  a  priori  separation  of  compact  and  diffuse 
layers  is  postulated  in  this  model,  and  it  also  does  not  involve  the  use  of 
local  dielectric  constants  (a  meaningless  quantity) .  The  calculations  show 
that  the  solvent  molecules  are  polarized  even  at  the  point  of  zero  charge. 

This  is  so  because  of  the  spillover  of  electronic  charge  into  the  solution, 
which  produces  non  zero  electric  fields  in  the  metal-solution  boundary  layer. 
Although  we  have  treated  this  effect  only  crudely,  the  recent  work  on  the 
GMSA  [16]  of  the  ion  dipole  mixture  allows  a  consistent  calculation  of  this 
effect. 


We  show  also  that  the  distance  dj  of  closest  approach  of  the  metal 
ions  and  the  solution  ions  is  a  crucial  parameter  in  the  determination  of 
the  differencial  capacitance  C.  However  this  is  also  a  drawback  since  the 
fact  that  dj  is  unknown  makes  the  comparison  of  C  with  experiment  a  not  very 
meaningful  excercise.  Indeed  a  value  of  dj  which  fits  the  experiments  could 
be  found. 


However,  the  model  that  we  have  discussed  here  represents  the  most 
comprehensive  one  that  we  are  aware  of.  We  are  undertaking  a  detailed  study 
of  the  influence  of  various  parameters,  such  as  ionic  concentration,  density 
and  dipole  moment  of  the  solvent,  etc.  on  the  differential  capacitance.  The 
results  will  be  published  in  the  near  future. 
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Table  I  -  Values  of  Che  profile  parameter  a  (Hg)  as  a 

function  of  d.  and  X  .  For  the  bare  surface  a  -  0.901  a.u. 
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16.6 

43.4 

Table  II  -  Differential  capacitance  of  the  interface  (uF/cm  ) 

as  a  function  of  dj  and  X.  The  concentration  is  c  =  0.1  M/1. 


FIGURE  CAPTIONS 


Fig.  I  :  Simplified  model  of  the  interface. 

6  position  of  the  ideally  smooth  wall. 
a  diameter  of  the  ions  and  solvent  molecules. 


Jellium  model  and  density  profile  of  the  dipoles  normal  to  the 
interface. 


Dependence  of  5^  on  the  parameter  dj  for  different  choices  of  A. 
(1)  Hg  A  -  15  -  (2)  Hg  A  -  30 

(3)  Hg  A  -  45  (4)  Ga  A  -  30. 


APPENDIX  1 


Electrostatic  Potential  Profile 

Integrating  (5)*  we  get,  (using  also  (I)  and  (3)  ) 

a)  For  negative  metal  charge  (z_  >  0) 

n  o 

V(z)  -  ea(2“*o>  (  -  as<  z  <  0) 

or 

V(z)  -  27m(z2  -  e+0t(z-zo)  f  o2  )  (0  <  z  <  z  ) 

o 

V(z)  -  2Tm[z  (2z  -  z  )  -  2/a2  ♦  (z  <  z  <  d.) 

O  O  O  1 

V(z)  -  2-rrn[z  (2d,  -  z  )  -  2/a2  +  e^^^/a2]  +  47m z JK 

o  i  o  o 

+ 

z  4ttuz  _ ,  ,  . 

+  4tt  /  dz’  P(z’) - °  e”ic(z-d^'  (z  >  d.) 

dl 

b)  For  positive  qu  (z  <  0) 

n  o 

V(z) - 2lTn  e^^^o)  (-  a»  <  2  <  z  ) 

a2  ° 

V(z)  =  -  27m  [(z  -  z  )2  +  2/a2  -  e-a(*-*0)/a2]  (z  S- 2  <  °> 

o  o 

V(z)  -  27rn  [z  (2z  -  z  )  -  2/a2  +  e~^(z~Zo)JaZ  ]  (0  <  z  <  d.) 

o  o  i 

For  d.  <  z  <  <»  the  equations  are  those  of  the  case  q._  <  0  . 
i  M 


APPENDIX  2 


We  give  explicit  formulas  for  the  coefficients  of  Eq.  23 

C.  -  —  1.6423  n2/3  +  0.25037  n1/3  +  0.006584] 

l  7m 
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Ihes£  coefficients  are  computed  once  and  for  all  for  a  given  metal  and  solvent 


